arXiv: 1507.05837vl [math.AP] 21 Jul 2015 


A contact problem for viscoelastic bodies with inertial effects 
and unilateral boundary constraints 


Riccardo Scala 

Dipartimento di Matematica, Universita di Pavia, 
Via Ferrata 1, 27100 Pavia, Italy 
E-mail: riccardo. scalaOunipv. it 

Giulio Schimperna 

Dipartimento di Matematica, Universita di Pavia, 
Via Ferrata 1, 27100 Pavia, Italy 
E-mail: giusch04@unipv. it 

July 22, 2015 


Abstract 

We consider a viscoelastic body occupying a smooth bounded domain H C R 3 under the 
effect of a volumic traction force g. The macroscopic displacement vector from the equilibrium 
configuration is denoted by u. Inertial effects are considered; hence the equation for u contains 
the second order term u tt- On a part T_d of the boundary of 12, the body is anchored to a support 
and no displacement may occur; on a second part T jv C <912, the body can move freely; on a third 
portion Tc C <91 2, the body is in adhesive contact with a solid support. The boundary forces 
acting on Tc due to the action of elastic stresses are responsible for delamination, i.e., progressive 
failure of adhesive bonds. This phenomenon is mathematically represented by a nonlinear ODE 
settled on To and describing the evolution of the delamination order parameter 2 . Following the 
lines of a new approach outlined in [5] and based on duality methods in Sobolev-Bochner spaces, 
we define a suitable concept of weak solution to the resulting PDE system. Correspondingly, we 
prove an existence result on finite time intervals of arbitrary length. 

Key words: second order parabolic equation, viscoelasticity, weak formulation, contact problem, 
adhesion, mixed boundary conditions, duality. 
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1 Introduction 

Let 12 C M 3 be a smooth bounded domain of boundary T. We assume 12 to be occupied, during 
the fixed reference interval (0,T), by a viscoelastic body. No restriction is assumed on the final time 
T > 0. The displacement vector with respect to the equilibrium configuration in 12 is noted by u. 
Hence, we may assume that u satisfies the following equation 

u tt - div(Ve(u t ) +Ee(u)) = g, in 12, (1.1) 

where, given a vector-valued function v, e(v) represents its symmetrized gradient; g is a volume 
force density; V and E are the viscosity and elasticity tensors, respectively, assumed symmetric, 
nondegenerate, bounded, and depending in a measurable way on the variable x £ 12 (we refer to 
Section [2] below for the precise assumptions). Inertial effects may occur and are mathematically 
represented by the second order term u tt . 

Let T := <912 be the boundary of 12 and let us assume T to be decomposed as T = r£>urArUrc, 
where the portions TV, for X = D,N, C, are assumed to be relatively open in T, mutually disjoint, and 
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to have strictly positive 2-dimensional measure. In addition to that, the distance between Tc and is 
assumed to be strictly positive. In other words, the contact surface Tc and the support Td need to be 
strictly separated by Tjv. In the portion Tp of the boundary, we suppose the body to be anchored to a 
rigid support; hence no displacement may occur, or, in other words, u satisfies a homogeneous Dirichlet 
condition. In the part T^v, we assume that the body is allowed to move freely, which corresponds to 
asking u to satisfy a homogeneous Neumann boundary condition (with standard adjustments we could 
equally consider the nonhomogeneous case, corresponding to the physical situation where a surface 
traction is exerted on Tjv). Finally, on Tc (here, C stands for “contact”), the body is assumed to be 
in adhesive contact with a hard surface, like for instance a wall. This configuration has two effects: 
firstly, the (trace of the) displacement uonTc must be directed towards the interior of f2. Namely, 
we ask the following constraint to be satisfied: 


u • n < 0 , on Tc, ( 1 . 2 ) 

where n is the outer normal unit vector to 12. The meaning of d is that the material in 12 may 
not penetrate the wall, but it may well detach from it. Secondly, the contact between the body in 12 
and the wall may cause “damage”, i.e., loss of adhesive properties. This phenomenon is often referred 
to as delamination process. Mathematically, it is described by means of a second (scalar) variable z, 
defined on Tc over the same time interval (0,T). The “bonding function” z takes the form of an order 
parameter ; namely, z(t 7 x) represents the fractional density of adhesive bonds that are active at the 
time t £ (0, T) and at the point x GTq- In particular, when z = 1, there is total adhesion, whereas 
when z = 0 the bonds are completely broken. A value z £ (0,1) denotes a partial loss of adhesivity. 
In addition to that, we assume that the bonds, once broken, cannot be repaired; namely, we require 
the time derivative z t to be nonpositive at each (t,x) £ (0, T) xTc. 

In view of the above description, only the values z £ [0,1] and z t £ (—oo, 0] are meaningful 
(or, as we will often say, are “physical”). Hence, the mathematical equation(s) for z must enforce 
in some way these constraints (or, equivalently, exclude the “unphysical” configurations). This scope 
may be reached by relying on the theory of maximal monotone operators (cf. the monographs [ 3 ) (HE 1 ). 
Indeed, the equation for z, settled on Tc, may take the form 

a(z t ) + z t +/3(z) B a - ^ |u| 2 , (1.3) 

where a and (3 are monotone graphs in Rxl entailing the required constraints. For instance, we may 
assume a = d/(_ 00 i o], i.e., the subdifferential of the indicator function of (—oo, 0 ], whereas we may 
take /3 = dl\o :+00 \, i.e., the subdifferential of the indicator function of [0, +oo). Hence, a provides 
the nonpositivity of zt (i.e., the irreversibility of damage), whereas /3 guarantees the nonnegativity of 
z (i.e., the fact that when the adhesive bonds are completely broken no further damage may occur). 
The constraint z < 1 (i.e., the fact that the material cannot be more than completely integer) is 
automatically guaranteed by zt < 0 once it is z < 1 at the initial time. 

The positive constant a on the right hand side of m ( more generally, we may also admit 
a to depend in a suitable way on i £ Tc) has the meaning of a threshold under which the elastic 
stresses are not strong enough to cause delamination. Namely, when A |u | 2 is less than a, the right 
hand side of El is nonnegative; hence no damage is created. Instead, delamination may occur for 
a — A |u | 2 < 0. Note also the term z t on the left hand side of & meaning we assume here the 
damage evolution to be rate-dependent. 

In order to get a closed system, one has to specify how the behavior of the damage variable 
influences the behavior of u. This is made explicit by the (up to now missing) boundary condition for 
u on Tc, which we assume in the form 

— (Ve(u t ) + Ee(u))n = y(u ■ n)n + zu. (1.4) 

Here, 7 is a third maximal monotone graph enforcing the non-penetration constraint (cf. Remark 12.II 
below for further comments). In particular, if we take 7 = a = c>/(_ 00 , 0 ], we obtain that (the trace 
of) u is directed towards the interior of 12, or, at most, we may have tangential displacements. In the 
sequel we will actually allow for more general assumptions on 7 (and this is why we decided to use a 
different symbol to denote it). 
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Models for contact, delamination and damage in elastic media are becoming very popular in 
the recent mathematical literature. The evolution law of the adhesive damage variable z and its link 
with the displacement on Tc is based on the so-called concept of Fremond delamination (see m, 
here considered with the presence of nonnegligible viscosity of the adhesive (whose effects arise from 
the term zt). Usually in the context of delamination models, the mechanical system consists of two 
(or more) elastic bodies attached upon an interface. Here we assume for simplicity to have only one 
body placed in Q and attached to the support in Tc; however we observe that our results could be 
extended to the general case by trivial generalizations. A dynamic model for delamination without 
viscosity of the adhesive is also considered in [3T] , where the evolution equation for the variable u is a 
variant of m which takes higher order stresses into account. A dynamic model where also thermal 
effects are considered has been analyzed in m- Other related models, among many, can be found in 
Enzmanaiiznisi- A model coupling (ED with an equation for z similar to ED including viscosity 
can be found in |22j . The main difference between our model and the one in [22l stands however in 
the presence of the unilateral constraint ED, which represents also the main mathematical difficulty 
occurring here. Actually, when one considers mechanical models with inertial effects (i.e., containing 
the second order term u tt ), enforcing the constraint ED by the methods of monotone operator theory 
usually gives rise to regularity issues; namely, interpreting monotone graphs as monotone operators 
in the usual L 2 (or L p ) framework is generally out of reach. For this reason, in most of the related 
literature, this difficulty is overcome by restating the equation containing the constraint as a variational 
inequality. This is the case, for instance, of the recent papers mm ■ In [2] , existence of a solution to 
a system closely related to (U.1D + (I1.4I) (hence, the damage variable is not explicitly considered there) 
is proved by discretization techniques. In m, well posedness of a system also accounting for the 
evolution of z is proved by restating both equations as variational inequalities. It is worth observing 
that, differently from here, in m the evolution of z is assumed in to be reversible, corresponding to 
the choice a = 0 in the equivalent of (11.31) : moreover, the quadratic term on the right hand side of 
the same relation is truncated (this simplifies the proof of existence compared to our case). 

In the present paper, we actually prefer to follow a partially different approach, based on 
our recent work [5] (in collaboration with E. Bonetti and E. Rocca), where a strongly damped wave 
equation for a real valued variable u containing a general constraint term is analyzed. The basic idea 
stands in writing a weak formulation where test-functions are chosen in a “parabolic” Sobolev-Bochner 
space V (in particular, in [5] , V = 7J 1 (0, T; L 2 (Q)) D L 2 (0, T; 7J 1 (fl))). In that setting, the monotone 
graph providing the constraint is restated as an operator acting in the duality between V and V'. 
Although it can be shown rather easily that this reformulation is factually equivalent to the statement 
as a (suitable) variational inequality, it presents a number of notable advantages: firstly, it clarifies 
the regularity of the constraint term, which usually corresponds to a physical quantity; secondly, it 
permits us to prove further properties of solutions, like the energy inequality, or to analyze the long¬ 
time behavior, which will be the object of a forthcoming work. Finally, in this setting we can still take 
advantage of the basic tools of monotone operator theory in Hilbert spaces in order to prove existence. 

It is worth pointing out a further mathematical difficulty of the present model compared to 
similar ones. Namely, the coupling term |u| 2 /2 in (11.31) has a sort of critical growth in space dimen¬ 
sion 3. This is due to the mixed boundary conditions complementing ED, where the “Neumann” 
part ED additionally accounts for the nonsmooth constraint term. In this setting, the best space 
regularity we can hope to obtain for u is H 1 . This translates into an irregularity of the trace on 
F and into a (maximal) irregularity of |u| 2 /2. In view of this fact, the identification of the terms 
a{zt) and j3(z) in the “doubly nonlinear” equation (11.31) is somehow nontrivial and requires a careful 
combination of monotonicity and semicontinuity tools. 

We finally point out that evolution equations with inertial terms and bilateral constraints, like 
for instance u ■ n = 0 on Tc, (corresponding to the case where only shear displacements are allowed) 
are mathematically much simpler to deal with. These models arise, for instance, in physical systems 
subjected to high pressure. Such a behavior is often referred to as “Mode II” evolution, in constrast 
to “Mode I” evolutions, that are those where the constraint is unilateral, as in ED- 

The remainder of the paper is organized as follows: in the next section we introduce our 
assumptions on coefficients and data and state a mathematically rigorous weak formulation together 
with our main result. This states existence of at least one weak solution of suitable regularity. The 
proof occupies the remainder of the paper. In particular, in Section [3] a regularized problem is 
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introduced and existence of a local solution to it is shown by means of a (Schauder) fixed point 
argument. Next, in Section [4j the approximation is removed by means of suitable a-priori estimates 
and compactness methods. Other qualitative properties of weak solutions, like the energy inequality, 
are also discussed there. 


2 Assumptions and main results 

We set H := L 2 (Q) and H := L 2 { fi) 3 . Moreover, for k > 1, we introduce 

H k D :={ueH k (n-,R 3 ):u = OonT D }, (2.1) 

and we denote by H\~ k the respective dual spaces. We put V := H l D , i.e., 

V := {u € H\Q;R 3 ) : u = 0 on Tq}. (2.2) 

We also set Vo : = and recall that V' Q = ff _ 1 (fi; R 3 ). The spaces V and Vo are seen as 

(closed) subspaces of /f 1 (f2; R 3 ) (and in particular they inherit its norm). The duality between V and 
V' will be indicated by (•, •). 

In the sequel we shall frequently use the continuity of the trace operator 

from H l (Vl) to H l ^ 2 {Tc) and from H 1 ^) to L 4 (Tq) (2-3) 

and its vector analogue. Moreover, the trace operator will be generally omitted in the notation; 
namely, functions defined in Q and their traces on Tq will be indicated by the same letters. 


Strong formulation. This can be stated as: 

u tt - div(Ve(u t ) + Ee(u)) = g, in(0,T)xO, (2.4) 

a{z t ) + z t +/3(z) 9 a - i|u| 2 , on (0,T)xT c , (2.5) 

- (V£(u f ) + Ee(u))n £ )(u • n)n + ;u, on (0,T)xTq, (2.6) 

complemented with the additional boundary conditions 

u = 0, on (0,T)xffl, (2.7) 

(Ve(u t )+Ee(u))n = 0, on(0,T)xTAf, (2.8) 

and with the Cauchy conditions 

u| t =o = u 0 , u t | t= o=ui, z\ t =o = zo, (2-9) 


where the first two relations are assumed a.e. in f2, while the third one is stated a.e. on Tq. 

Remark 2.1. It is worth noting that we considered the homogeneous condition (12.81) just for the 
sake of simplicity. Indeed, one could deal with the case of a nonzero boundary traction f on the 
right hand side with standard modifications. Moreover, it may be also worth emphasizing that, if 
<r := Ve(u t ) + Ee(u) denotes the stress tensor, then the boundary relation (12.61) could be split into 
its normal and tangential parts as follows: 

— fi n := — (ern) • n G y(u • n) + zu n, on (0, T ) X Tq, (2-10) 

— er t := — ern + ((ern) • n)n = ^u t := z( u — (u • n)n), on (0,T) X Tq. (2.11) 

This corresponds exactly to the conditions considered, e.g., in 0 (1.23-24)] for v = 0, i.e., when no 
friction is assumed to occur on Tq. 

As explained in the introduction, we are not able to provide a solution to the strong (pointwise) 
formulation of the system. Consequently, we need a weaker notion of solution. To introduce it, we 
start with stating our assumptions on coefficients and data: 
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(a) The tensors V, E € L°°(f2;R 81 ) satisfy, a.e. in O, the standard symmetry properties 

Y ijkl — Y jikl = Y ijlk ~ Y klij 5 E ijkl = E jiki = Ejy/fc = E klij • (2.12) 

Moreover, V, E are assumed to be (uniformly in f 1) strongly positive definite; namely, there 
exists a constant k > 0 such that 


E(x)e : e > k\e\ 2 , 

for a.e. x £ f 1 and any symmetric matrix e £ 
inequalities, for any v £ V there holds 


V(x)e : e > k\e\ 2 , (2-13) 

R 9 . Hence, in view of Poincare’s and Korn’s 


/ Ee(u) : e ( v ) > k||u||^, / Ve(u) : e(u) > k||u||^, (2-14) 

if! Jo. 

for a (possibly different) constant k > 0. 

(b) We set a = cU(_oo.o]) moreover, we assume /? and 7 be maximal monotone graphs inRxM such 
that D(/ 3 ) = [0, +00) and Dfy) = (— 00, 0]. In particular, an admissible choice is /3 = <9/; 0 ,+oo) 
and 7 = oo.o]- We recall that the domain D(b ) of a graph b C K. x R is the set {r £ 
R : b{r) ^ 0 }. We denote as 3 , / 3 , 7 suitable convex and lower semicontinuous functions from 
R to (—00, +00] such that a = da , /3 = d/ 3 , 7 = cfy. Hence, we have in particular a = /(_oo,o]. 
We also suppose that 


/3(r),7(r) > 0 resp. for all r £ D((3),D( 7). 


(2.15) 


(c) The initial data satisfy 

u 0 £ V, 7 (u 0 • n) £ L 1 (rc), ui £ H, (2.16) 

together with 

£0 G T°°(r c ), 0 < zo < 1 a.e. on r c , /3°fyo) € L 2 (Y G ). (2.17) 

Here, for r £ D(/3), /3°(r ) is the element of minimum absolute value in the set /3(r) (cf. [9]). 
Using the definition of subdifferential one may easily prove that (12.171) implies in particular 

/3fyo) £ L\T c ). (2.18) 


(d) The volumic force satisfies g £ L 2 (0, T; V'). 

(e) We let H be a smooth and bounded domain of R 3 . We assume that T = dfl satisfies T = T G U 
TjvUTc, where Tjc, for X = D, N, C, are relatively open in T, mutually disjoint, and have strictly 
positive 2-dimensional Hausdorff measure. Moreover we assume the distance d(Tc,r.D) > 0 
and we require that Tc is smooth as a subset of T and has at most finitely many connected 
components. Specifically, we assume that the boundary of Tc in T is an at most finite union of 
curves of class C 1 . 


Remark 2.2. It would also be possible to consider more general choices for a. For instance, we 
may ask a to be a maximal monotone graph satisfying D(a) = (—00, 0], the analogue of (12.151) . plus 
some additional conditions regarding the behavior near 0. However, this would give rise to a number 
of technical complications in the proof, whence we decided to restrict ourselves to the basic choice 
a = <9/(_ oc , 0 ] ■ On the other hand, our somehow general assumptions on /3 and 7 do not require any 
additional technical work. 

Remark 2.3. Assumption (12.151) essentially states that /? and 7 must have some coercivity at 00. For 
what concerns /3 this is in fact just a technical assumption, in view of the fact that we will prove that 
z < 1 almost everywhere. However, it may be useful in the approximation. Note that the analogue 
of (12.151) also holds for a = /(_oo,o]- 
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Energy functional. System (12.4I) - (I2.6D has a natural variational formulation. Namely, it can be 
seen as a generalized gradient flow problem for a suitable energy functional. It is worth pointing out 
this structure from the very beginning. To this aim, we will obtain the energy estimate directly from 
the system equations. Of course, such a procedure has just a formal character at this level since we 
have not yet specified which is our notion of solution and the related regularity. That said, we first 
test (I2.dll by u t . Using also (12.1211 with the boundary conditions (12.711 on T^ and (12.811 on Tat, it is 
easy to obtain 


^ L (^ Ut l 2 + i!> IEf ^ U ) 1 £( ' U 0 + X V£< ' U *' ) :£( ' U *' ) = / ((Ve(u t )+Ee(u))nj •u t +(g,u t ). (2.19) 
Then, we test m by Zt and integrate over Tc- A simple integration by parts in time gives 

<^ / ^ + + / ( a ( Zt ) + z *) Zt = f 2 ( u ' u *)- ( 2 - 20 ) 

Next, to cancel the terms on the right hand sides of (12.19l) - (12.20l) . one scalarly multiplies (12.61) by — u f 
and integrates, to obtain 


'To 


((Ve(u t ) + Ee(u))n) 


u t = - 


s/ r u (un) -/ rc 2(uu *>' 


Hence, taking the sum of (12.1911 . (12.2011 . (12.2111 . we (formally) obtain the energy identity 

d 


cl t 


£(u, u f , z) + V(ut,z t ) = (g, u t ), 


with the energy functional £ = £(u, u t , z) given by 
I I ^|u t | 2 + ^Ee(u) : e(u)'j + j 


T c 


/3(z) - az + -z\u\ 2 +y(u ■ n) 


and the dissipation integral(s) V = T>(u t ,Zt) dehned as 


V ■= 


:= / Ve(u t ) : e(u t ) + / (a(z t ) + z t )z t . 
h 2 Jv c 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


(2.24) 


Note that, in view of assumptions (a)-(d), both £ and V enjoy suitable coercivity properties. Observe 
also that the right hand side of (12.2211 accounts for the contribution of external volumic forces. 

Although relation (12.2211 corresponds to a basic physical property of the model, the mathemat¬ 
ical procedure we used is formal under many aspects. The main problem is related to the occurrence 
of the nonsmooth multivalued graphs a, /?, 7. Hence, we will see in Thm. [7UT1 below that, for weak so¬ 
lutions, we will be only able to (rigorously) prove a weak version of (12.2211 in the form of an inequality 
(cf. (12.4411 below). 

In order to state our precise concept of weak solution we start with introducing some more 
functional spaces: 


V := (2.25) 

U~H\Q,T-,H*{T C )), (2.26) 


and we let V' and %' be the respective dual spaces. Moreover, for all t € (0,T], we set 

V t :=H\Q,t;V), (2.27) 

(2.28) 

with the dual spaces V' t and T-L' t , respectively. Note that H is exactly the space of traces (on Tc) 
of the elements of V (and similarly for Ht and V t ). In the sequel, we shall note as ((•, •)) the duality 
pairings with respect to both space and time variables. For instance, that symbol may note the duality 
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between V and V' or also that between H and . When working on subintervals (0, t), t < T, we 
will use the notation ((•, -}} t (e.g., that may denote the duality between V* and V' t ). Not to weight up 
formulas, we will use the symbol (•, •) for the scalar product in both H and L 2 (Tc)- The norms in 
H and L 2 (Vc) will be sometimes simply noted by || ■ ||. The double brackets ((•, •)) will represent the 
L 2 -scalar product in time-space variables (for instance, in L 2 (0,T; H) or in L 2 (0, T; L 2 (Tc))). On 
time subintervals, we will use the notation ((•, •)) t . 

Next, we define the convex functional 

G : L 2 (0,T-, L 2 (Tc)) —t [0,+oo], G(v) := f f 7 ( 1 ;). (2.29) 

J 0 JYc 

Then, if one considers the subdifferential dG in the (Hilbert) space L 2 (0,T; L 2 (Tc)), it is well known 
that this coincides with the realization of the graph 7 . Namely, for v,t) £ i 2 (0,T; L 2 (Tc)), one has 

•q £ dG[v) q(t,x) £ "/(v(t,x)) a.e. on (0,T) x Tc. (2.30) 

Hence, in particular, v complies with the constraint represented by 7 . On the other hand, in our 
specific situation, we will not be able to interpret 7 in the above sense, due to regularity lack coming 
from the occurrence of the term u tt . For this reason, following the lines, e.g., of 0 (cf. also w , we 
provide a suitable relaxation of 7 . To this aim, we identify L 2 (0, T; L 2 (Tc)) with its dual by means 
of the natural scalar product, obtaining the chain of continuous and dense inclusions 

n c i 2 (o,T;i 2 (r c ))~i 2 (o,r;i 2 (rc))'cH'. (2.31) 

Then, the above constructed Hilbert triplet permits us to relax 7 in the following sense: we define as 
7 W (where the subscript “re” stands for “weak”) the subdifferential of the restriction of G to the space 
T~L with respect to the duality between H and TL'. Namely, for v £H and 77 £ 71', we set 

q £ lw{v) (( 7 , w — v)) + G(v) < G(w) for all w £ T~L. (2.32) 

A precise characterization of 7 ^, which is a maximal monotone operator from T~L to 2 W , is carried out 
in 0 Sec. 2 ] following the lines of results first proved in .ID] (see also 0 )- Here we just mention the 
fact that, for any v £ H, there holds the inclusion y(v) C 7 m(v), which may however be strict. On the 
other hand, once we know that 77 £ 7 ^( 77 ), then v is still necessarily almost everywhere nonpositive 
(hence, it satisfies the constraint); moreover, the inclusion 77 £ j w (v) has a precise “measure-theoretic” 
interpretation in terms of the original graph 7 (see 0 Sec. 2] for more details). The analogue of 7 ^, 
noted with the same symbol for the sake of simplicity, may be constructed also on the space %t, i.e., 
working on time subintervals. 

We are now ready to introduce our concept of weak solution: 

Definition 2.4. Let T > 0 and let Assumptions (a)-(e) hold. We say that (u, 77 , z, £ 1 , £ 2 ) is a weak 
solution to system m-m ^ 


u G W llOO (0, T; H) n H l { 0, T; V), 

(2.33a) 

u t £H\0,T- V')nW(0,T;iT^ 2 ), 

(2.33b) 

z £ W 1 ’ oo (0,T; L 2 (T c )), 

(2.33c) 


(2.33d) 

6,6 £L°°(0,T;L 2 (T c )) 

(2.33e) 


and the following properties are satisfied: 

(i) For all tp £ V it holds 

(u t (T),<p(T)) - ((u + ((Ee(u), £(<£>))) + ((¥e(u t ), e(ip))) 

+ «¥ F ' n» + f [ zu-ip = ((g,ip)) + (u 1 ,tp(0)), (2.34) 

J 0 J Tc 
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with the initial conditions (EHD. Correspondingly, for every t £ [0,T) there exists r £ TLt such 
that 


(u t (i), ¥>(*)) - ((u t ,¥> t )) t + ((Ee(u), e(ip))) t + ((Ve(u t ), e(<p))) t 

+ ((m , ■ n)) t + [ [ zu-<p = ((<g,tp)) t + (u 1 ,<p(0)), (2.35) 

J 0 JYc 

for all tp £ V t . Moreover, the functionals r) and rj^ are compatible, namely, if p £ Vt satisfies 
p(t) = 0 a.e. on f l, then we have 


{{V(t),<P ' n)) t = ((v, F ■ n)), (2.36) 

where Tp represents the trivial extension of tp to V (i.e., Tp(s,x) = 0 for s £ [t,T\ and x £ ft), 
(ii) For a.e. t £ (0,T) there holds 

6(0 +zt{t) +£,i(t) = a- ^|u(t)| 2 , a.e.onT c . (2.37) 

(in) The following graphs inclusions hold true: 

&G0Ca.e. on (0,T) xr c , (2.38) 

6 6 a(^t), a.e. on (0,T)xr c , (2.39) 

£ 7w(u • n). (2.40) 

Moreover, for all t £ (0, T) we have 

1(t)6«((u'n)L ( o,i)). (2.41) 

Note that, since z and the term with rj in (12.3411 are concentrated on Tc, if we choose <p £ i7 1 (0, T; Vo) 
in (12.3411 and integrate by parts in time, we get back 

(u tt ,p)+ [ Ee(u):e(y)+ f Ve(u t ) : e(<p) = (g,ip), (2.42) 

J Q JQ 


a.e. in t £ [0,T], where the first duality product makes sense in view of the first (I2.33bl) . Hence, in 
particular (12.3411 implies (II.Ill in the sense of distributions. 

Let us now recall that the energy of the system was defined in (12.231) . Moreover, analogously 
with (12.241) . we introduce the energy dissipation as 

D := [ Ve(u t ) : e(u t ) + [ (6 + z t )z t . (2.43) 

J n Jr c 

We are now ready to state our existence theorem, constituting the main result of the present paper: 

Theorem 2.5. Let T > 0 and let Assumptions (a)-(e) hold. Then there exists at least one weak 
solution (u, 77 , z, 6 , 6 ) t° Problem (12.4D - (I2.9I) . in the sense of Def. 12.41 Moreover, for all times 
t 2 £ [0, T] and for a.e. ti £ [0, ^ 2)7 the following energy inequality holds: 

£(h)+ f 2 T>{-) < £(ti) + [\g, u t ). (2.44) 

J t\ Jt 1 


The proof of the above result will occupy the remainder of the paper. 














3 Regularized problem and local existence 

In this section we introduce a regularized version of our problem and prove existence of a local (in 
time) solution by means of a suitable fixed point argument. In view of the fact that this procedure is 
standard under many aspects, we omit most details and just present the basic highlights. 

First of all, for e £ (0,1) intended to go to 0 in the limit, we take suitable regularizations a E , 
/3 E and of the maximal monotone graphs a, /3 and 7 . In particular, we will consider the Yosida 
approximations defined, e.g., in [9], to which we refer the reader for details. Here we just recall that a e , 
j3 e and y E are monotone and Lipschitz continuous functions defined on the whole real line. Moreover, 
they tend, respectively, to a, /? and 7 in a suitable way, usually referred to as graph convergence in 
1x1 Namely, for any [r; s] £ a and any e £ (0,1), there exists r e £ R such that r e —>• r and 
a E (r E ) — > s in K as e \ 0, with analogous properties holding for /3 and 7 . In view of our choice 


a = dl(- oc, 0 ], we may explicitly compute 

a e (r) = £ -1 (r) + , A e := (Id+a £ ) _1 (r) = -( r)~ H-—r( r ) + - (3-1) 

e + 1 

Notice that A e is a nonexpansive operator. 

Our aim will be to solve, at least locally in time, a regularized statement, which, in the strong 
form, can be written as follows: 

u tt -div(Ve:(u t ) + Ee(u)) = g, in H, (3.2) 

a e (z t ) + z t + /3 e (z) = a - ^|u| 2 , on r c , (3.3) 

- (Ve(u t ) + Ee(u))n = 7 E (u • n)n + (z) + u, on Tc, (3.4) 


coupled with the initial conditions (12.911 and the boundary conditions (I2.7M2.8I1 . The function ( z) + in 
(13.41) denotes the positive part of z. Actually, z is not guaranteed to be nonnegative at the approximate 
level since the (smooth) function /3 E does not enforce any constraint. 

In fact, we will deal with a weak formulation of the above system, to which we will apply 
Schauder’s fixed point theorem. Hence, we start with introducing the fixed point space: for a small 
but otherwise arbitrary number s £ (0,1/2), and for T 0 £ (0, T] to be chosen at the end, we set 

S s (T 0 ) := {u £ H s (0,T 0 -,H) nL 2 ( 0 ,T o ;F^ s (fi;K 3 )) : v| r £ T 4 ( 0 , T 0 ; T 4 (r ; R 3 ))}. (3.5) 

The space S s (To) is naturally endowed with the graph norm, noted as || • ||$,(t 0 ) for brevity, which 
turns it into a Banach space. Notice also that the trace of v makes sense in view of the assumed 
U 2 + s -space regularity. Then, for some M > 0, we consider the closed ball 

B m := {v £S s (T 0 ) : |M|s s (t 0 ) < M}. (3.6) 

Note that the choice of M > 0 is essentially arbitrary. Its value will in fact influence the resulting 
final time To, but this is irrelevant at the light of the subsequent uniform estimates. 

The basic steps of our fixed point argument are carried out in the next three lemmas. 

Lemma 3.1. Let u 0 satisfy (12.1611 and zq satisfy (12.1711 . More precisely, let us set 


U := ||u 0 ||v + ||ui||jj + || 7 (u 0 • n)|| £ i(r c ), (3.7) 

z '■= ll-olk^rc) + ll^(^o)IUi(r 0 )- (3-8) 

Let also u £ Bm- Then there exists one and only one function z, with 

IMU-(o,T 0 ;.L 2 (r c )) + IMIrri(o,r 0 ;L 2 (r c )) < Q(e - \ Z, M), (3.9) 

satisfying, a.e. on (0, To) xTc, the equation 

(z t ) + z t + (3 s {z) = a — - |u| 2 , (3.10) 


with the boundary condition z\t=o = zq- Here and below, Q denotes a computable nonnegative-valued 
function, increasingly monotone in each of its arguments, whose expression may vary on occurrence. 
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PROOF. Using (l2.17D - fl2.18H with the graph convergence /3 s /3, it is not difficult to prove that, 

at least for e £ ( 0 , 1 ) small enough, there holds 

||zo|U 2 (r 0 ) + ll^ e (^o)||Li(r c ) < 2Z. (3-11) 

Then, existence of a solution z can be proved by using standard existence results for ODE’s. The 
regularity (13.91) can be inferred simply by testing (13.101) by z t . Actually, integrating over Tc, we then 
obtain 

Ir P E ^ + Jr aE ^ Zt + = J r 

< ^ll-t|li 2 (r c ) + C ( 1 + ll n lll 4 (r c ;R 3 ))- (3-12) 

Observe that the integration by parts is allowed in view of the smoothness of (3 e . Then, by the 
first m, the second integral on the left hand side is nonegative. Hence, using also (13.111) . we infer 

ll-tlU 2 (o,r 0 ;L 2 (rc)) + ll^ £ ( z )IU o °(o,T 0 ;i 1 (r c )) <Q{M,Z,s 2 ), (3.13) 

whence, using again (13.111) to estimate z from Zt, we obtain (13.91) . Finally, to ensure uniqueness of z 
one can use standard contractive arguments. For instance, if z\ and Z 2 are two solutions, one may test 
the difference of the corresponding equations (13.101) by the difference {z\ — z^)t and use monotonicity 
and Lipschitz continuity of of and /3 e together with Gronwall’s lemma. We omit details. I 


Lemma 3.2. Let Uo, zq, M, u as above. Let also z be the function provided by the previous lemma. 
Then there exists one and only one function u, with 

ll u llw 1 .°“(0 1 To;.H') + IMIifqo.TojV) + ll u IU oo (0,T 0 ;V) < Q{ £ \ Af, Z, t/), (3-14) 

satisfying, a.e. on (0, To) and for any p £ V, the equation 

(u tt,<p) + [ (Ve(u t ) + Ee(u)) : e(p) + [ ( 7 E (u • n)n + (z) + u) • p = (g, <p), (3.15) 

Jn J r c 

together with the Cauchy conditions u|t=o = uq and u t |t = o = ui. 


PROOF. The weak formulation (13.151) is obtained from (13.21) simply testing by cp £ V, integrating by 
parts, and using the boundary condition m- Then, existence of at least one solution can be proved 
by adapting the procedure given, e.g., in [7]. Here we just reproduce the corresponding regularity 
estimate, which is needed in order to get (13.141) . Namely, we take p = u t in (13.151) . and, proceeding 
as in the Energy estimate detailed before, we get 


d_ 

df 


£ (^Iwl 2 + ^Ee(u) :e(u)) + ^ q e (u ■ n) + ^ Ye(u t ) : e(u t ) 

= - / (-) + ( u ' u t) + ( 5 ,u t ) < ||-||L2 (ro) ||u|| L 4 (rc) ||u t || i 4 (rc) + ||ff||v'||ut||v 

JT C 

< Qie- 1 , M, Z)\\u\\ v \\u t \\v + ||t/||v'||ut||v < Q(e _1 , M, Z)||u||^ + c||g||^ + ^||u t ||^, 


(3.16) 


where k is the same constant as in (12.131) . Note that, to deduce the last inequalities, we used (13.91) 
together with the trace theorem (cf. (12.31) 1 and Young’s inequality. Now, as before, from (13.71) and the 
graph convergence y e —> 7 , we can easily prove, at least for e £ ( 0 , 1 ) small enough, 


||u 0 ||v + ||ui||j* + ||7 e (u 0 • n )||ii ( r c) < 2 U. (3.17) 

Hence, integrating (13.161) in time, recalling (12.13D - fl2.14D . and using Gronwall’s lemma, we readily 
obtain that u satisfies (I3.14D . 

To obtain uniqueness, we can proceed similarly as above. Namely, we assume to have two 
solutions Ui and U 2 to (13.15D . take the difference of the corresponding equations, and substitute 
p = (ui — 112)4 therein. Then, the Lipschitz continuity of y E and the properties of the trace operator 
permit us to get a contractive estimate via a procedure similar to (13.16D . I 
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Lemma 3.3. Let uo, zq, M, u as above. Let also z be the function provided by Lemma 1 1 and u 
be the corresponding function provided by Lemma 13.21 Let us consider the map 

T : B m —t kF 1,oo (0, To; H) fl T°°(0, To; V), T : u u. (3.18) 

Then, at least for e £ (0,1) sufficiently small, we can take T 0 £ (0, T\, possibly depending on e, M, U 
and Z, such that the map T 

(a) takes values into Bm; 

(b) is continuous with respect to the (strong) topology of S s (To); 

(c) maps Bm into a compact subset of <5^ (To). 

PROOF. (a) Thanks to (13.141) . to (12.31) . and to standard embedding and trace theorems, we have 

ll u llL 4 (0,To;L 4 (r;R3)) < cT^ ‘ Ju|| L oo( 0! T 0 ;L 4 (r;R3)) < cT 0 1 / 4 ||u|| i o 0 (o i T 0 ;V') — C ^C(3.19) 
ll u lliJ s (0,T 0 ;.H') < c||u||j f i(o i T 0 ;H') < c Tq^ 2 || u ||vF 1 - <x, (0,T 0 ;ff) < cTq^Q, (3.20) 

H U llL2(0,T 0 ;ff3 + s (O;R3)) - C ll U IU 2 (0,T o ; V) < cTp 7- || u|| L = (0 ,T o; v) < cT^Q, (3.21) 

where we wrote Q in place of Q(e~ 1 ,M, Z, U ), for brevity, and where c > 0 are embedding constants 
independent of Tq. Hence, we can choose To sufficiently small, possibly depending on e, so that 

IMk(To) < cT^Q(£-\M,Z,U) < M, (3.22) 


as desired. 

(b) Let {u,,} C Bm and let u„ — >■ u in S s (Tq). Let also z n and z be the corresponding functions given 
by Lemma [Til let u„ = T(u„) and let u = T(u) be the corresponding solutions given by Lemma lT2l 
We have to prove that 

u n —^ u strongly in <S s (To). (3.23) 

First, repeating, with the proper adaptations, the uniqueness argument sketched in Lemma 13.11 (and 
using in particular the Lipschitz continuity of a e and /3 s ), we can easily show that 

lim || z n - zWhi^ t . L 2(y p = 0. (3.24) 

n/*oo 

Next, we work on equation (13.151) . Proceeding similarly with (13.161) and performing standard manip¬ 
ulations, it is not difficult to obtain 

lim (||u„ — u|| W i,oo(o,T 0 ;/f) + ll u n ~ u ll L“ (0,T o ; V) ) = 0, (3.25) 

n/'oo 

This relation, also on account of (12.31) . implies (13.23L as desired. 

(c) The proof is similar to the above one, but a bit more tricky. Indeed, we still consider a sequence 
{u„} C Bm, but we now just assume that 

u„ -»• u weakly in S s (To). (3.26) 

Then, with the same notation as above, we need to show that at least a subsequence of {u„} satisfies 
(13.231) . To prove this fact, we first observe that, thanks to standard interpolation and (compact) 
embedding results, there exists a (non-relabclled) subsequence of n such that, for some p £ (1,2) 
depending on the choice of s, 

u„-s>u strongly in T 2p (0, T 0 ;T 2p (r c ;K 3 )), (3.27) 

whence we have in particular 

|u„| 2 -»• |u| 2 strongly in L p (0, T 0 ;L P (T C )) and weakly in T 2 (0, T 0 ; T 2 (r c ))- (3.28) 
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Let us now write (13.1011 for the index n and for the limit (where z is the solution corresponding to u), 
and take the difference. Rearranging term and applying the inverse operator (Id+a £ ) _1 , we get the 
relation 

(z n - z)t = (Id +a e )~ 1 ^ — i|u„| 2 + i|u| 2 - f3 e (z n ) +/3 e {z)Y (3.29) 

Then, testing by \z n — 2 | p_1 sign(z n — z) and using the Lipschitz continuity of (Id+ctd) -1 and of /T 
together with Holder’s and Young’s inequalities, it is not difficult to arrive at 

J)^ Zn ~ - C JW Z ™ “ z \\lp(Tc) + H n « - n llL2p (rc )l|u n + u||P 2p(rc) ). (3.30) 

Hence, applying Gronwall’s lemma, we obtain 

z n ^z strongly in L p (0,T o -L p (T c )). (3.31) 

Moreover, by the Lipschitz continuity of (Id+cY) -1 and of /3 £ , we can easily prove that the analogue 
of (13.311) holds also for ( z n ) t . Repeating the a priori estimates given at point (b) to obtain (13.2411 . we 
then conclude that 

z n —} z weakly in H^O, T 0 ; T 2 (Tc))- (3.32) 

Combining (13.311) and (13.321) we also get 

(z n ) + (z) + strongly in 47(0,T 0 ;L r (r c )) for every r £ [1,2). (3.33) 

Next, considering equation (13.151) with right hand side depending on z n and repeating the usual energy 
estimate, it is easy to obtain 

u„ —» u weakly star in W 1,oo (0, To; H) D IL 1 (0, T; V) (3.34) 

for some limit function u. In particular, due to standard compact embedding results for vector-valued 
Sobolev spaces, this entails 

u„-s-u strongly in iL s (0, T 0 ; iT) CL 2 (0, T 0 ; H5 +S (H;R 3 )). (3.35) 

Hence, to get (13.231) . it remains to show that 

u„ —> u strongly in T 4 (0,T o ; T 4 (T;K 3 )). (3.36) 

This will be proved at the end. Preliminary, we may notice that u solves the limit equation (13.151) . 
Actually, we have 

(z r j) + u„ —>• (z) + u weakly in 17(0, T 0 ; L r (Tc\ R 3 )) for some r > 1, (3.37) 

as a consequence of (13.331) . (|3.35l) . and of the continuity of the trace operator from H^ +S (Q) to H S (T) 
for any s £ (0,1/2). Hence, it turns out that u = T(u). To conclude the proof, we then need to show 
(13.361) . More precisely, we will reinforce (13.341) proving that 

u„ —> u strongly in LL 1 ’ oo (0, T 0 ; H) D iL 4 (0, T; V). (3.38) 

Then, (13.361) will follow from (12.31) . To get (13.38[) we need to use a semicontinuity argument, which we 
just sketch. Indeed, the same procedure will be repeated in the next section under more restrictive 
assumptions. We actually test (13.21) . written for u„, by u„, and integrate. Then, using and 

performing some integration by parts, we arrive at 

f / Ee(W) : £(u e ) + \ [ ¥e(u £ (t)) : e(u £ (f)) = \ f ¥e(u 0 ) : e(u 0 ) - f f 7 e (u £ • n)u £ ■ n 
Jo Jn z Jn z Jn Jo Jn 

+ [ [ l u t| 2 - «(*),u e (t)) + (ui,u 0 ) - [ [ (z e ) + \u £ \ 2 + f {g, u £ ). (3.39) 

Jo Jn Jo Jt c Jo 
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We also test 02D. written for the limit u, by u, obtaining an analogue relation. Then, we take the 
limsup of (13.3911 at the level n and we compare the outcome with (13.3911 written for u. Treating the 
terms on the right hand side by owing to the smoothness of y £ and by semicontinuity tools (see the 
next section for details), we then obtain, for every t £ (0,T o ], 

limsup f f f Ee(u E ) : e(u E ) + ^ f Ye(u £ (t)) : e(u E (f)) 

e\0 \io if! 2 Jq 

< f f Ee(u) : e(u) + i f Ve(u (t)) : e(u(*)). (3.40) 

J 0 J Q ^ J £2 

Thanks to the symmetry and coercivity properties of the tensors E and V (cf. assumption (a)), we 
then get (13.38H . whence (13.3611 . Summarizing, we have 

T(u„)h>T(u) strongly in S s (T 0 ), (3.41) 

which actually holds for the whole sequence u„. Hence, the map T is compact. This concludes the 
proof of the lemma. I 

As a consequence of the three lemmas, we can apply Schauder’s fixed point theorem to the map T, 
at least for £ £ (0,1) sufficiently small. This provides existence of at least one local solution (u E , z £ ) 
to the approximate system. To be precise, ( u £ ,z £ ) satisfies, a.e. in (0, To), and for any cp £ V, 


( u tt^) + [ E ^(u e ) :e(¥>)+ / Ve(uf) : e{p) + [ q E (u E • n)<p ■ n + [ (z e )+u e • p = (g, p). (3.42) 
Jn J n J r c Jt c 

Moreover, (13.311 holds a.e. on (0, To) and the couple (u E ,z E ) also complies with the initial condi¬ 
tions (12.91) . 

Remark 3.4. It is worth remarking that (u E , z e ) also satisfies an approximate version of the energy 
equality. Indeed, comparing terms in (13.4211 . one can easily prove that u E £ H 2 (t),Ty V'). In partic¬ 
ular, this is sufficient in order for p = uf £ T 2 (0,To; V) to be an admissible test function in (13.4211 . 
Hence, testing also (13.31) by and proceeding as in the last part of Section [21 we may infer 



uf,z E )+T E (u E ,* E ) 


{di u t} ~ f |(- e ) u E -u E , a.e. in (0,To). 
JTc ^ 


(3.43) 


where £ e and D £ are the approximate energy and dissipation functionals and the last term appears 
in view of the occurrence of the positive part in Note that here all integrations by parts are 

fully justified in view of the fact that (3 and 7 have been replaced by their regularized counterparts. 
Integrating (13.4311 in time, we then also obtain the additional regularity 


u c £C' 1 ([0 ) To];JT). 


(3.44) 


4 Global existence for the original system 

In order to prove Theorem 12.51 we will show that, as e —> 0, the regularized solutions (u e , 2 : e ) con¬ 
structed before tend, in a suitable way and up to the extraction of a subsequence, to a weak solution 
to the original problem. It is worth noting from the very beginning that, at least in principle, the 
functions (u E ,z e ) are defined only on some subinterval (0,To) possibly smaller than (0,T) and also 
possibly depending on e. However, in view of the fact that we shall derive a set of a-priori estimates 
that are independent of To, standard extension arguments imply that (u e , z e ) can in fact be extended 
to the whole of (0,T), and the same will hold for the limit solution ( 11 , 2 ). Hence, in order to reduce 
technical complications, we shall directly assume with no loss of generality that (u £ , 2 e ) are defined 
over (0, T) already. In particular, the approximate energy equality (13.4311 and the related regularity 
(13.44[) turn out to hold on the whole of (0, T). 

That said, we proceed with the proof, which is subdivided into various steps, presented as 
separate lemmas. 
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Lemma 4.1 (Extension of boundary functions). Let xp £ 14 1 / 2 (rc; R 3 )- Then, there exists p =: 
Rxp £ V such that <p\r c = xp in the sense of traces. Moreover, the operator R : xp 1-4 p is linear and 
continuous from H 1 J 2 (Tc', R 3 ) to V, namely there exists cn > 0, depending only on 12, Tc, Tjv, T^j, 
such that 

IMIv < c IIV’llff 1 / 2 (r c ;R 3 )- (4-1) 

PROOF. We first observe that, in view of the regularity assumptions (e) on 12, using extension by 
reflection and cutoff arguments, xp can be extended to a function xp defined on the whole of T in such 
a way that xp is 0 a.e. on To and 

(r ; R 3 ) < c\\M\hv 2 (r c ;R 3 )- (4.2) 

Moreover, we can build the extension in such a way that the map xp i—>■ xp is linear. As a second step, 
we construct p as the solution to the elliptic problem 

— Atp = 0 in 12, <p = xp on T. (4-3) 

Then, the desired properties follow from (14.21) and standard elliptic regularity results. I 


Lemma 4.2 (Step 1: first a priori estimate). There exists a constant M > 0 independent of e £ (0, 1) 
and of T 0 £ (0,T] such that 


Moreover, we have 


and, for all t £ (0, T), 


uE ||w 1 .' x, ( 0,T;JT) < Af, 

(4.4a) 

ue ||fl’ 1 (0,T;V) < M, 

(4.4b) 

u ttlU 2 (0,T;V^) < 

(4.4c) 

u t bf(o,T;R d 2 ) — 4/, 

(4.4d) 

Ur 1 (o,T;L 2 (r c )) < M. 

(4.4e) 

u " ■ n)||Li(o,T;Li(r c )) < M , 

(4.4f) 

U e • n)\\ w < M, 

(4.4g) 

|| 7 e (ii e • n)\\ H[ < M. 

(4.4h) 


Proof. Let us note as /? E and suitable antiderivatives of (3 e and y E , respectively. Then, using 
the standard relation |6 E (r)| < |5°(r)| holding for any maximal monotone graph b C Rxl (cf. [1]), and 
recalling assumption (12.151) . it is not difficult to prove that, up to suitable choices of the integration 
constants, we can assume /3 E and to be nonnegative. 

That said, let us take p = uf in (13.421) . multiply (13.31) by d t (z e ) + (i.e., the time derivative 
of the positive part of z e ), sum the resulting expressions, and integrate in time. Noting as H the 
Heaviside function, we then observe that 

[ j3%z E )d t (z E ) + = f p e (z°)H(z s )d t (z s ) = 1 f B e (z e ), (4.5) 

Jr c Jrc dt Jr c 

where the function B e (r) coincides with /3 E (r) for r > 0 and is identically equal to /3 E (0) for r < 0. 
Hence, we infer 

^IKWII 2 + \ [ Ee(u E (t)) : e(u E (t)) + f f ¥e(u E (t)) : £ (uf(t)) + [ 7 E (u E (i) • n) 
z z in Jo Jn Jv c 

+ \ [ (zT(t)\u E (t)\ 2 + f \\d t (z e )+\\ 2 + f [ a e (Zf)d t (z e ) + + f B%z%t)) 

^JT C J 0 J 0 JT C 

= a [ [ d t {z e ) + + [ (g,u E ) + i||ui|| 2 + i / Ee(u 0 ) : e(u 0 ) + / 7 e (u 0 • n) + f B E (z 0 ) 

Jo Jr c Jo z z Jn Jr c Jr c 

+ \J r toluol 2 <c+i^ \\dt(z s ) + \\ 2 + ^ || Vu E || 2 +c||fif|| 2 2 (0iT;V , ) , (4.6) 
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where, in the last line, we have used Young’s, Poincare’s and Korn’s inequalities together with the 
analogue of (Id.111) . (Id. 171) . and the fact that 


1 

2 



2o|u 0 | 2 < ||zo||i 2 (r c )|l u o|li 4 (r c ;R 3 ) < c||2o||z, 2 (r c ) ll u o ||v < c - 


Hence, recalling (12.1dl) . we easily obtain (14.4al) and (14.4bl) . 

Let us now test ES by zf. Then, proceeding as before and noting that 


(4.7) 




— \\ Z t IU 2 ( 0 ,T;L 2 (r c )) ll u ~ lll 4 ( 0 ,r;L 4 (r c ;R 3 )) 

— 9 \\ z t\\ 2 L 2 (0,T-,L 2 (r c )) + C II U ^ lli 4 ( 0 ,T;V)> 


(4.8) 


using (I4.4bl) to estimate the last term, we readily infer EH . 

Next, we choose p £ Vo in (Id. 421) . Then, the boundary integrals go away and a simple 
comparison of terms permits us to obtain (I4.4cl) . 

Let us now prove ( |4.4g[ ). To this aim, let ip £ T-L. Then, in view of the regularity assump¬ 
tions (e), we can think the outer unit normal n to H on Tc to be the trace of a smooth function, denoted 
with the same symbol, defined on an open neighbourhood A C i 3 of Tc. As a consequence, ip := ipn 
lies in iL 1 (0, T; iL 1 / 2 (Tc; R 3 )). Hence, applying Lemma l4~Tl we may construct p = Rip £ iL 1 (0,T; V) 
such that p\v c = V 1 ) an d 


||‘Pllff 1 (o,T ; v) < cll^llffqo.itffV^rcsR 3 )) S c||V’||-h- (4.9) 

Taking such a p in (ld.42L integrating in time, performing suitable integrations by parts, comparing 
terms, and applying (I2.1dl) and Korn’s inequality, we then have 


|(( 7 e (ir -n),V>))| < ||uf|| L 2 (0iT:fJ) || ¥ 5 t || i 2 (0jT;ir) + ||u|(r)|| J i||vj(r)|| JJ + ||ui||jj||v3(0)||h 

+ c ll Vu t llL 2 (0,T;ff) ||V<p||i 2 (Q iT;f/ ) + c||Vm ||i2(0,T;ff)l|V¥ , llL2(0 > T ; ff) 

+ ll(~ £ ) + u " llL 2 ( 0 ,T;L 2 (rc;R 3 ))ll < Pl|i 2 ( 0 ,T;L 2 (r c ;M 3 )) + I|9||l 2 (0,T;V') I|¥ > IIl 2 ( 0 ,T;V) 

< c||^||v < cIIV'II-h, (4.10) 


where the constants are provided by (l4.4aD - d4.4cD . (|4.4ejl . and the continuity of the trace operator. 
Hence, ( |4.4g[ ) follows. Repeating the same argument on subintervals, we also obtain d4.4hD . 

Finally, to get d4.4fl) let us apply Lemma 1X11 to the function ip = n. As noted above, we may 
assume n £ L7 1 (0, T; iL 1 / 2 (Tc; R 3 )). Hence, we obtain an extension m £ iL 1 (0, T;V) of n to the 
whole domain fb Let us now plug p = u £ + m into dd.42D . Then, estimating the resulting right hand 
side as in d4.10D . it is not difficult to arrive at 


7 e (u £ • n)(u e + m) ■ n < c||u e + m||v < c, 


(4.11) 


0 JTc 


the last inequality following from d4.4bD and Lemma 14.11 applied to n. Then, we claim that the left 
hand side can be estimated from below as follows: 

[ [ 7 e (u e -n)(u £ + m)-n = [ [ 7 E (u £ -n)(u £ -n + 1) > K 0 ||7 e (u e •n)|| z ,i (0iT:i i (ro)) -c. (4.12) 
JO JTc: Jo JVn 


with kq > 0 and c > 0 independent of e. To prove the last inequality in d4.12D we can use the 
assumption D( 7 ) = (— 00 ,0] with the coercivity d2.15D . Indeed, on the one hand this implies that, for 
r > — 1 / 2 , there holds 

7 e (r)(r + l) > ^| 7 e (r)|. (4.13) 

On the other hand, due to d2.15D . either linv^-oo y(r) < 0 (and the same holds for 7 E ), whence we 
can reason as in d4.13D also for r << 0, or it is linv^-oo 7 (r) = 0 (and, again, the same holds for j e ), 
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so that there is nothing to prove because in that case the L 1 -norm of 7 e (u e ■ n) may only explode on 
the set where u e • n > —1/2. Hence, we have proved (Id. 121) . which, combined with (14. lit) , gives (I4.4fl) . 

Finally, we take ip £ H 2 D in (13.421) . Then, noting that ||v 5 llL“( ro ) — cllvllff 2 ,; a comparison of 
terms in (13.421) . together with estimate (I4.4fl) . permits us to obtain (I4.4dl) . which concludes the proof 
of the lemma. I 

Lemma 4.3 (Step 2: second a priori estimate). There exists a constant M > 0 independent of 
£ £ (0,1) and of To £ (0,T] such that 

ll/3 e (z E )l|L“(0,T;L2(r c )) < M, (4.14a) 

ll a£ (-t)llL=°(o,T;L 2 (r c )) <M, (4.14b) 

ll~t llL“(0,T;L 2 (r c )) < M. (4.14c) 

Proof. Let us multiply equation (13.101) by -^/3 E (z e ) = (^ £ )'(z e )zf, so to obtain 

— / r(^)i 2 + / m\z e M\ 2 + f m^)a E {zt)zt 

^Jrc jTc Jrc 

= ^ t J r (a-\\u E \ 2 )f{z E ) + {vf-ulWiz*), (4.15) 

whence, integrating over (0,f), 0 < t <T, and using that |/3 e (-)| < |/3°(•)| together with assumption 
(12.171) . we get 

\\\P E (z E ml H r c ) +J*J r G8 e )V)l*t | 2 + J* l m’(z E )a%z E t )zt 
<lW(zo)\\l Hrc) + 

+ ll u iU“(0,T;L 4 (r c ;R 3 )) (2 H U f lli 2 (0,T;L 4 ( rc ;R3)) + lll 2 (0,T;L 2 (r c ))) 


< c(l + \\P S (z e )\\ 2 L2(o ,T;L 2 (r c )))’ (4-16) 

where we have also used Young’s inequality and the estimates (14.41) . Now, since (/3 s )' > 0 and a e 
is monotone and satisfies a E (0) = 0, we may use Gronwall’s lemma to obtain (I4.14al) . Moreover, 
applying the nonexpansive operator (Id+cY)^ 1 , we may rewrite (13.101) in the form 

z t = (Id+a £ ) _1 (a — ^|u £ | 2 — P e {z E )^J. (4.17) 


Comparing terms in (14.171) , we then get (|4.14cl) . With this information at disposal, we go back to (13.101) 
and a further comparison argument gives also (I4.14bl) . I 

Lemma 4.4 (Step 3: converging subsequence). There exist limit functions (u, z, 77 , £ 1 , £ 2 ) such that, 
for a (non relabelled) subsequence of e —>• 0, there holds 

u E —^ u weakly in i? 1 (0,T;V) and weakly star in kF 1,oo (0,T; H), (4.18a) 

uf — *• u t weakly in L7 1 (0,T;V'o) and weakly star in BV(0,T; Hfj 2 ), (4.18b) 

z € — y z weakly star in W 1,oo (0,T; L 2 (Tc)), (4.18c) 

7 e (u e ■ n) —^ 77 weakly in hi! , (4.18d) 

f3 e (z £ )^fi weakly star in L°°(0,T;L 2 (Tc)), (4.18e) 

a e (zt) —>■ £ 2 weakly star in L°°(0, T; L 2 (Tc)), (4.18f) 

together with 

U( —» u t strongly in L 2 (0,T;H), (4-19) 

Uj(t) — 1 u t (t) weakly in H , for all t £ [0,T]. (4.20) 
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Moreover for all t £ (0, T) there exists r?( t ) £ TL' t such that, for the same subsequence of e \ 0 
considered before, 


7 e (u e -n) L (o weakly in U' t , (4.21) 

PROOF. Convergences (I4.18all . (|4.18bl) . (I4.18dl) follow from estimates (14.41) . Convergences (I4.18el) 
and (I4.18H) follow from (14.141) . Moreover, (I4.14cl) implies (I4.18cl) . Let us now show (14.191) and (14.201) . 
Thanks to (14.18al) and (I4.18bl) . we can apply the generalized form of Aubin-Lions lemma ((221 Corollary 
4], [23 Corollary 7.9]) with the triple of spaces V CC H C Vq. This provides (14.191) . To see (14.201) . 
we first observe that such weak limit holds true in the space V' 0 by (I4.18bl) . Then the claim follows 
thanks to the fact that ||u t (i)||jj < M for every t £ [0,T] by (14.18al) and (13.441) . 

It remains to show (14.211) . Let us set v e := ( z e ) + . Then, thanks to (14.18cl) . there exists a 
nonnegative function v such that 

v e —^ v weakly in L 2 (0,T; L 2 (Tc))- (4.22) 

Moreover, by (14.18all - (|4.18b|l . the Aubin-Lions lemma, and (12.31) . we infer 

u £ —y u strongly in 7/(0, T;//(Tc; R 3 )), for all r £ [1,4). (4-23) 

Combining (I4.18al) . (14.221) and (14.231) . we obtain 

w £ u £ vu weakly in L 4/3 (0, T; L 4/3 (T C ; R 3 )). (4.24) 

Now, for all t £ (0, T), let if £ Tit, let if := if n and let ip be the extension of if provided by Lemma ITTl 
Then let us define a functional r/( t ) £ as follows: 

dV(t),ip))t ■= [ (ut^t) - (u t {t),cp(t)) + (ui,y>(0)) - [ [ Ee(u) : e(<p) 

Jo Jo Jn 

Ve(u t )-.e(ip)-f [ vu ■ <p + f (g,<p). (4.25) 

J o J r^7 J o 

Now, thanks to (14.18al) and (14.221) . it is seen that the value of ((?/(*), V , ))t is exactly the limit of 

[ f 7 e (u £ • n)if = f (u e tl ip t ) - (uf(t),ip(t)) + (u £ ,ip(0)) - f f Ee(u e ) : e(<p) 

Jo Jr c Jo Jo Jn 

u £ u e • ¥>+ [ (g,<p), (4.26) 

Jo 

which is obtained integrating (13.421) in time, rearranging terms, and performing some integration by 
parts. Hence, in particular is independent of the extension map if e-t ip. We also observe that, 
thanks to (I4.18al) - (|4.18bl) and (14.241) . the right hand side of (14.261) converges (to the right hand side 
of (14.251) 1 with no need of extracting further subsequences. The thesis follows. I 


f [ Ve(uf) : e(<p) - f f 
J o J J o *^r c 



Lemma 4.5 (Step 4: refined convergence for u). There hold the additional strong convergences: 

u E —y u strongly in L 2 (0,T;V), (4.27) 

u £ (t) —> u(t) strongly in V for all t £ [0, T]. (4.28) 


Moreover, the functions r) and r /( t ) are identified as follows: 

V G 7 m(u ■ n), 

V(t) e 7i«(( u ' n ) L (o,t))- 


(4.29) 

(4.30) 
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PROOF. Let us define, analogously with (12.2911 . the convex functional 


G s : L\ 0, T; L 2 {Y C )) -t [0, +oo), G e (v) := f [ (4.31) 

J 0 J T c 

Then, as observed before, the subdifferential dG e in L 2 (0, T; L 2 (Tc)) factually coincides with the 
graph 7 £ . On the other hand, we may interpret the function y E also as a monotone operator from 
% into W. Indeed, if v £ H, then it results 7 E (i>) £ L 2 (0,T; L 2 (Vc)) C H' thanks to the Lipschitz 
continuity of y e . Moreover, for all it, v £ H we have 

((7 e («) - l e {v),u - u)) = (( 7 e (u) - 7 e {v),u- v )) > 0. (4.32) 

Also, for u £ 77, 7 e (w) belongs to the subdifferential 7 ^ of (the restriction to "H of) G E at the point u 
computed with respect to the duality pairing between H 1 and 77. Indeed, we have 

7 E (u)(w - u) < [ [ (j e (v) - 7 e (u)) = G e (v) - G E (u), (4.33) 

JO JTc 

for all v £ 77. In short terms, we have y E C 7 ^,, where 7 ^ acts as a maximal monotone operator 
between 77 and 77'. 

Now, thanks to the monotonicity with respect to e of the functionals G E , owing to [2 Theorem 
3.20] and [2 Theorem 3.66], the maximal monotone operators 7 ^, converge to 7 ^ in the graph sense: 

V [x; y\ £ 7 ™, 3 [x e \y e ] £ 7 ^ such that [x e ; y e ] -» [x; y], (4.34) 

where the convergence is intended with respect to the strong topology of 77 x 77'. 

In order to take the limit of the equation for u, we first observe that (13.4211 . after integration 
in time over (0,T), can be equivalently rewritten in the form 

((Ee(u E ), e(<*>))) + ((Ve(u E ), e(<p))) + (( 7 e (u e -n),<p- n)) 

= ((u E , ¥ 3 t ))-(u E (T), V (T)) + (u 1 , V (0))- [ T [ v £ u £ .p + ((g,p)), (4.35) 

J 0 J Yc 

for any function p £ V, where v £ = {z s ) + . Let us also notice that the counterpart of (14.351) on 
subintervals (0 ,t) could be stated analogously. We now aim to let e \ 0. To start with, we take 
p = u E in (14.351) . Then, integrating by parts and rearranging terms, we get 

((7 e (u e ■ n),u E ■ n)) = -((Ee(u E ),e(u E ))) - ^ f Ye(u E (T)) : e(u E (T)) + ]- j V£(u 0 ) : e(u 0 ) 

J Q J 

+ (( u tX)) - K( T )X( T )) + (ui,u 0 ) - [ f u £ |u E | 2 +((g,u E )). (4.36) 

J 0 J Tq 

Now, recalling (14.231) - (14.241) . we may take the limit e \ 0 in relation (14.351) . Actually, using also 
(14.18al) and (I4.18dl) . we easily arrive at 

((Ee(u),e(v>))) + ((Ve(u t ), e(<p))) + ((rj,p- n)) 

= (K,<X) - (u t (T),^(T)) + (ui,<p(0)) - f f vu-p+ ((g,<p)). (4.37) 

J 0 JTc 

In particular, in view of the fact that u £ V, we may take p = u. Proceeding as for (14.361) . we then 
have 

((? 7 ,u • n)) = — ((Ee(u), e(u))) - i / Ve(u(T)) : e(u(T)) + i f Ve(u 0 ) : e(u 0 ) 

JTl J Q 

T 

+ ((u t ,Ut)) - (u t (T),u(T)) + (ui,u 0 ) - [ [ v\u\ 2 + ((g } u)). (4.38) 

Jo JTc 


((7 e (u),v-u))= f ( 

Jo Jr 
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Now, we use the so-called Minty’s semicontinuity trick in order to identify the function 77 . To this 
aim, we take the lirnsup as e \ 0 in (14.3611 and we compare the outcome to (14.3811 . Here, the key 
point stands in dealing with the integral term. Actually, setting 

$:lxl 3 4 l, $0,u) := 0) + |u| 2 , (4.39) 

we may observe that $ is continuous and nonnegative. Moreover, it is convex in u for any v £ R. 
Hence, Ioffe’s semicontinuity theorem (cf., e.g., [14]) yields 

$(u, u) < liminf f f $(u £ , u e ) = liminf / [ u e |u E | 2 . (4.40) 

Jo Jr c e N0 Jo Jr c 

Indeed, we already know that (u) + = v and (u E ) + = V s . 

With (14.401) at disposal, taking the lirnsup of (14.361) . using (I4.18al) . (|4. 191) - (14.201) . (14.231) and 
semicontinuity of norms with respect to weak convergence, and finally comparing with (14.3811 . we get 

lirnsup (( 7 e (u e • n), u E • n)) < (( 77 , u • n)). (4-41) 

£\(0 

Thanks to the fact that 7 E (u E • n) £ 7 ^,(u e • n) and to the graph convergence (14.341) . this suffices to 
prove (14.291) . i.e., to identify 77 (cf. [TJ Prop. 3.59], see also [9J Prop. 2.5]). 

At this point, we are able to reinforce the strong convergence of u E in (14.231) . To this aim, let 
us first notice that, in view of (14.341) . there exists a sequence [x e \y e ] £ y E such that [ x e ; y e ] —> [u n; 77 ] 
strongly in TL x TL'. Then, by monotonicity of 7 ^, we have 

« 7 e (u E ■ n),u E • n)> > « 7 e (u £ • n),x e )) + «t, e , u E ■ n)) - ((y £ ,x £ )). (4.42) 

Then, taking the liminf, noting that all terms on the right hand side pass to the limit, and using 
(14.41 [) . we may conclude that 

lim il e (u E ■ n), u e ■ n)) = (( 7 ?, u ■ n)). (4.43) 

erNjO 

At this point, we rearrange terms in (14.361) rewriting it in the following way: 

((E£(u e ),£(u e ))) + i / V£(u e (T)) :£(u e (T)) = -(( 7 E (u e -n),u E -n)) + i / V£(u 0 ) : e(u 0 ) 

J Q 1/ 

+ ((u E ,u E ))-(uf(T),u E (T)) + (u 1 ,u 0 )- [ T [ u £ |u E | 2 + ((g, u E )). (4.44) 

J 0 J Tq 

We also rearrange terms in (14.381) in a similar way. Then, taking the lim sup of ()4.44[) , using (14.431) , 
and comparing with (the rearranged) (14.381) . we obtain 

lirnsup ("((E£(u E ),£(u E ))) + i [ Ve(u e (T)) : e(u s (T)) \ < ((Ee(u), e(u))) + ]- f Ve(u(T)) : e(u (T)). 
e\,o \ ^ Jn J J Jn 

(4.45) 

Using assumption (a), we then get (14.271) . Then, repeating the argument on subintervals (0,t), we 
also obtain (14.281) . As a further consequence, combining (14.281) with the uniform boundedness (I4.4bl) . 
we get 

u e —x u strongly in I/ 4 (0,T; V"), whence in i 4 (0, T; L 4 (rc; K 3 )), (4.46) 

the latter property following from (12.31) . Relation (14.46[) will play a key role when we take the limit 
of the equation for z. 

To conclude the proof of the lemma, we need to show (12.411) . This, however, follows simply 
by repeating the above argument on the subintervals (0 ,t), for t <T. I 

Lemma 4.6 (Step 5: refined convergence for z). There hold the additional strong convergences 

z £ — y z strongly in L 2 (0,T-L 2 (Tc)), (4.47) 

z e {t) —>• z{t) strongly in L 2 (Tc), for all t £ [ 0 ,T], (4.48) 
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PROOF. We prove (14.471) and (14.481) by adapting an argument due to Blanchard, Damlamian and 
Ghidouche [5[ Lemma 3.3] (see also [12]). To start with, applying A e (cf. (13.11) '). equation (13.101) may 
be equivalently rewritten as 

zl = A s (a-\\u s \ 2 -F{z*j). (4.49) 

For S € (0, i), we subtract from (14.491) the analogue expression with S in place of £, then we multiply 
the result by 2 r e — z 5 . After integration on (0, t) x Tc, we get 


\\\z E (t)-z s (t)\\l Hrc) 

= / / (M a - 

J 0 Jt c v 


-F( z t))-A e (a- l -\ 

u E 2 — /3 s (z s ))^(z E 

-z s ) 

+ // (M a - 

J 0 JY C y 

i|W| 2 -P\z s ))- 

A e (a - 

||u 5 | 2 -0V)))(* e 

-z*) 


+ J' J (-W- 

Jo Jv c v 

l\u s \ 2 -p s (z s ))~ 

As (a - 

V| u5 | 2 -f3 s (z s )))(z e 

-z s ) 


=: I\(t) + 12 (f) + I 3 (t) 





(4.50) 

By the nonexpansivity of 

A e , we first compute 




0 

^_ 

VI 

K "1 

|u e + u 5 ! u E — u 5 | 

z e -z s \ 




< ||u E - u 

L<5 lL 2 ( 0 ,t;L 4 (r c ;R 3 )) 

u E + u 5 

L“(0,(;L 4 (r c ;R 3 )) “ 

^ <5 ||L 2 (0,t;L 2 (r c )) 


< c||u e - 

,i<51| 2 

u llL 2 (tM;L 4 (r c ;R 3 )) 

+ lk e - 

z ^\ L 2 ( 0 ,t;L 2 (r c ))' 


(4.51) 

Next, recalling (13.11). 







m TT - JTi) (“ - - W-4) V - A 


10 JT C 


5+1 

,5II2 


< c(e + <5)(l + ||u , '||^oo(o it;i 4( rc;R 3)) + ||/3' 5 (^ 5 )||i=°(o,t;i 2 (rc)))ll- e _ A'|Ui(o,t ; L 2 (r c )) 

< c(e + <5)||z E - ^||i,i(o,t ; L 2 ( ra )), (4.52) 

the last inequality following from the previous uniform estimates. 

In order to estimate I\[t) we argue as in [T21 p- 270]. Namely, we use the identities 

z e -£p e (z E ) = R e (z e ), (4.53) 

P e (z e ) € P(R s (z E )), (4.54) 

valid for all e G (0,1), where R e is the resolvent of (3 at step e (see [U p. 27-28]). We have 


hf) = I 1,1 + h,2 ■= 


A e (a - ±|W| 2 - pfz E )) - A e (a - ±|u E | 2 - f3 s (zf) 


o Jr c ,„ pe{zt)-p 5 (z s ) 

x (R e (z E ) - R s (z s ))(p E (z E ) - p\z s )) 


+ /V (A e (a-|u E | 2 -^ E ))-A £ (a-i|u1 2 -^(^)))(^ £ (^)-^V)), (4-55) 

j o j r c 

where Tc.o is the set where the denominator does not vanish. Now, Zip is easily seen to be nonpositive 
thanks to the monotonicity of A e and to (14.531) . The second term I 12 is controlled as follows: 

1^1,21 < \\/3 £ {z e ) - P 5 (z s )\\ L 2 i0 ^L 2 {rc ))(\\£P e (z s ) - Sp s {z s )\\ L 2 ^ t . L 2 (rc)) ) <cf + S), (4.56) 

where we used (14.14al) twice. Gathering the estimates (14. 5 111 - (14. 561) . and setting <j>f) := || z £ (t) — 
z 5 (t) || i 2 ( ro ) we obtain, for every t € (0,T] and every s £ (0,f], the differential inequality 


<^ 2 (s) < ^-(^e + 5 + ||u" - u 5 || 2 2 ( 0 t;i 4( r(7 . R 3)) + ||z e ^ z' 5 ||j, 2 ( 0jt;i 2 ( rc p^ + c 2 (e + 5) J </>(•), (4.57) 
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for suitable ci, C 2 independent of e and 5. Then, applying the generalized Gronwall lemma Lemme 
A.5], we deduce 

1 2 

4>(s) < ci (e + S + ||u e - u ' 5 |li2( 0it;i 4 (rc;R 3)) + ||z e — ^ 5 |lL 2 (o,t;L 2 (r c ))) + c z( £ + ^) s > (4.58) 

whence, squaring, and choosing s = t, 

< C 3 + S + ||iT - u' 5 ||| 2 ( 0 it;i 4( rc;R 3)) + J (4.59) 

Applying once more Gronwall’s lemma, now in its standard form, and recalling (I4.46[) . we then arrive 
at 

||z £ (t) — z s (t) ||l 2 (p c ) —> 0 as |e + < 5 | —> 0, (4.60) 

for every t £ [0, T], which implies (14.4711 and (14.4811 in view of (I4.18cl) . I 

Lemma 4.7 (Step 6: limit flow rule). The limit functions provided by Lemma 14.41 satisfy condition 
(ii) of Def. 12.41 Moreover the inclusions (|2.38D and (12.3911 hold. 


PROOF. Using (|4.18cfh (I4.18cl) . (I4.18fl) . and (14.461) . we can take the limit in equation (13.31) and get 
back (12.3711 . or, in other words, condition (ii). Hence, it remains to identify and £ 2 - Firstly, we 
observe that inclusion (12.381) follows by combining (14.4711 with (14.1 Sell and using, e.g., [4j Prop. 1.1, 
p. 42]. Next, to prove (12.3911 . we use once more Minty’s trick (cf., e.g., (H Prop. 2.5]); namely, we need 
to check that 


lirnsup / / a E {zl)zl < 

£—>0 J 0 JY c 



(4.61) 


Then, we multiply (13. Kill by zf to obtain 

£ l ae(4K = -|l<llL 2 ( 0 ,T;Z, 2 (r c) )-^ ?(z e {T))+£ r(zo)+^ T ^ (a-i|u e | 2 )z?. (4.62) 


Then, taking the lirnsup, we may observe that the four terms on the right hand side can be managed, 
respectively, by (14.18cl) and semicontinuity, by (14.18cl) with the Mosco-convergence of (3 e to jd (cf. [T] 
Chap. 3]), by the monotone convergence theorem, and by combining (I4.18cll with (14.461) (note that 
having strong convergence of u £ in L 4 is essential at this step). As a consequence, we then infer 


lirnsup/ 1 [ a £ (zt)zt <-\\z t \\l H0T . L2iTc)) - [ P{z(T)) + f f}(z 0 ) + [ [ (a-^|u| 2 )^, 

J 0 JTc JT c «/ 0 jYc ^ 

(4.63) 

Now, testing (12.3711 by z t and applying the chain rule formula [9] Lemme 3.3, p. 73] to integrate the 
product term £1 z t , we see that the above right hand side is equal to f Q T f r f-iZt- Hence, (14.611) follows. 
In turn, this implies (12.391) . which concludes the proof. I 


Lemma 4.8 (Step 7: condition (i)). The functions u, z, 77 , and obtained in Lemma |4.4I satisfy 
condition (i) of Def. 12.41 


Proof. The statement follows from all the results obtained so far. In particular, we notice that 
the function v in (14.2211 coincides with (z) + , thanks to the strong convergence (14.4711 . Moreover, as 
a consequence of (12.38[) . we have z > 0 a.e. on (0,T) x Pc. This implies that ( z) + = z. Hence, 
we may take the limit as e \ 0 in equation (14.351) . which gives back (12.341) . Analogously, repeating 
the procedure on a subinterval (0 ,t), we obtain (12.351) . Finally, taking a test function <p £ Vt with 
ip(t) = 0 and considering its trivial extension Tp to the whole of (0, T), plugging Tp in (14.3511 and <p in 
the analogue of (14.3511 over the interval (0,t), and finally letting e \ 0, we obtain (12.3611 . I 

Lemma 4.9 (Step 8: energy inequality). The energy inequality (12.4411 holds for almost every t\ £ [0,T] 
and for every t^ £ 
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PROOF. We know that the approximate solution (u e , z e ) satisfies the energy equality (13.4311 . Inte¬ 
grating it in the time interval [ti,t 2 ], we get 

S e (t 2 )+ f V e (-) = £ E (h) + [ (g, u() - / f (z £ )-\i s -u £ t , (4.64) 

Jt\ Jt\ j t\ c 

for every 0 < t\ < t 2 < T. Now, we take the liminf on both hand sides. Then, standard semicontinuity 
arguments and the convergence relations proved so far permit us to prove that 

s s (t 2 )+ /‘Vo 

■It! 

for every ti, t 2 . Moreover, in view of the fact that z > 0 almost everywhere, it is easy to check that 
( z e )~ tends to 0 strongly in L 2 (0, T; L 2 (Fc)), whence the last integral in (14.641) vanishes in the limit. 
Moreover, it is also easily seen that the term depending on g passes to the limit. 

Hence, it just remains to control £ £ (t i), which is the more delicate point. Indeed, what we 
want to do is taking the lim inf of 

£ e (ti) = / (||u|(ii)l 2 + ^(u £ (H)) : e{u e (h)j) + [ ( - az £ {t 1 ) + ±z £ (h)\u £ 

J £2 j r q 

+ [ (^(^(H))+7 £ (u £ (ti)-n)). (4.66) 

JT C 

Then, in view of (14.27l) - (14.28l) . (14.191) . (14.461) and (14.18cl) . the liminf of the first two integrals is in fact 
a true limit and is given by the same couple of integrals rewritten without the e. On the other hand, 
this holds just for almost every ti; indeed, to take the limit of the integral of |uf (ii)|, we need to use 
(14.191) : in other words, we deduce a.e. (in time) convergence from L 2 -convergence by extraction of a 
subsequence. 

Finally, we need to control the last integral in (14.661) . This can be done by following closely 
the argument devised in [5J Sec. 3], to which we refer the reader. Actually, by that method we may 
prove more precisely that 

limsup / (/3 £ (z e (ti)) + 7 £ (u £ (fi) • n)) < J (^(z(h)) + 7 (u(*i) • n)), (4.67) 

again, for almost every t\. This concludes the proof. I 

Remark 4.10. In addition to (12.44|) . we may notice that the energy functional defined in (12.231) is 
lower semicontinuous with respect to the variable t € [0,T]. Indeed, as a consequence of (12.33aD - 
(12.33cl) . we have 


(4.65) 


rt 2 


£(t 2 ) + j £>(•) < liminf ^ 


ueC([0,T];n u t eC w ([0,T\-,H), 2 eC([0,I];L 2 (T c )). (4.68) 

Hence, the lower semicontinuity of t i —> £ (u(t), z(t)) follows easily from the convexity of /3 and 

7 and from lower semicontinuity of norms with respect to weak convergence. 
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